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We study steady-state properties of inelastic gases in two-dimensions in the presence of an 
energy source. We generalize previous hydrodynamic treatments to situations where high and low 
density regions coexist. The theoretical predictions compare well with numerical simulations in the 
nearly elastic limit. It is also seen that the system can achieve a nonequilibrium steady-state with 
asymmetric velocity distributions, and we discuss the conditions under which such situations occur. 
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I. INTRODUCTION 



Granular materials such as sand and powders have gen- 
erated much interest of late. Such an ensemble of par- 
ticles with macroscopic size is challenging since it may 
behave as a solid, a liquid or a gas. Size separation, pat- 
tern formation, avalanches, compaction and convection 
are just a few examples of the wide array of observed 
phenomena jlj. 

Flow underlies most of these phenomena and therefore, 
theoretical studies so far focus on formulatinga hydro- 
dynamical description appropriate to sand [|2|j2|]. These 
theories, stemming from the Boltzmann Equation, de- 
pend on the assumption of "molecular chaos", i.e., the 
assumption that no interparticle correlations exist. This 
assumption is far from obvious. As a dissipative dynam- 
ical system, a granular system has attractors in its phase 
space, which may cause correlations between particles. 
Under certain conditions, these attractors lead to a sin- 
gularity, inelastic collapse, which can not be explained by 
hydrodynamics. In one dimension, these attractors are so 
strong that hydrodynamics breaks down for the entire pa- 
rameter space — in a confined geometry, all particles, 
save one, form a practically stationary clump against the 
elastic wall, while the remaining particle moves rapidly 
back and forth between the clump and the heated wall. 
Such a state clearly violates partition of energy. 

In this study, we investigate the corresponding situa- 
tion in two-dimensions We consider inelastic hard 
spheres in a box where one wall is kept at a fixed tem- 
perature and the other three are reflecting (see Fig. [I]). 
Energy input at the heated wall balances the dissipation 
due to inter-particle collisions and the system can achieve 
a steady state. Unlike the one-dimensional situation, the 
density and the temperature profiles are smooth func- 
tions of the distance from the heated wall. In the steady- 
state, the momentum balance equations imply that the 
pressure is constant throughout the system. Particles 
move faster close to the energy source and more slowly 
deeper inside the system due to energy loss. Thus the 



density is greater farther away from the wall to maintain 
a constant pressure. 

Density variations may be a consequence of the above 
mechanism or an effect of the intrinsic attractors in 
the system which build up correlations among parti- 
cles. Assuming that the particles are not coherent in 
the quasi-elastic limit, we derive a differential equation 
that describes the density variation throughout the sys- 
tem. Heuristic expressions for the mean-free path, the 
equation of state, and the thermal conductivity are incor- 
porated into the energy flux balance equation to obtain a 
closed equation describing the steady-state density pro- 
files. The theory compares well with simulational results 
over a wide range of densities in the quasi-elastic limit . 
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FIG. 1. The system of interest. The length of the system 
is L, the width is W. The heated wall is at x = 0; the elastic 
walls are at x — L, y = 0, and y = W. 



Our theory is appropriate for a steady-state which is 
very close to an equilibrium state. We expect this to be 
applicable only in the nearly elastic limit. We argue that 
in the complementary situation of stronger inelasticity, 
hydrodynamics may still be relevant. To analyze the na- 
ture of the steady-state, we study velocity distributions 
throughout the system. We observe that the velocity dis- 
tributions exhibit scaling. This observation is then used 
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to obtain a qualitative description of the behavior of the 
system. 



II. THE QUASI-ELASTIC LIMIT 

Granular materials are different from ordinary fluids 
or gases in that the diameter of a particle may be com- 
parable to the collisional mean free path. To formulate 
a hydrodynamic theory, it is necessary to describe how 
the mean free path, the pressure, and the thermal diffu- 
sivity depend on the grain diameter, the number density, 
and the temperature. In this section, we obtain heuristic 
expressions for these quantities. 

Let us denote the number density of grains by p, the 
mean free path by I, and the grain diameter by d. With- 
out loss of generality, we set the particle mass and the 
Boltzmann constant to unity: m = ks = 1- The gran- 
ular temperature can be defined as the average kinetic 
energy per particle: T = ^(v 2 ). This quantity is well 
defined for an equilibrium state in which the particles 
have a symmetric velocity distribution. In such a case, 
v = \/2T, can be used as an approximate value for the 
average grain speed. In our system, quantities such as 
temperature, density, etc., are position dependent, and 
there is not a global equilibrium. However, in the quasi- 
elastic limit, we expect the system to be very close to 
local thermal equilibrium. 

Let x be the distance from the heated wall. In the 
steady state, all quantities vary only in this transverse 
direction. The energy balance equation can therefore be 
written 



dq 

dx 



0, 



(1) 



where q is the energy flux in the x direction. The sink 
term / accounts for the energy lost per unit area per 
unit time due to inelastic collisions. The energy flux is 
induced by a temperature gradient, q = —ndT/dx, where 
the coefficient of thermal diffusivity, k, is positive. Con- 
sequently, we find 

Tx( K dx-)= L (2) 



A. The Energy Sink 

Since collisions between grains are inelastic, kinetic en- 
ergy is continually transferred into heat. For simplicity 
we neglect rotation and thus, the degree of inelasticity 
can be parameterized by < r < 1. When two par- 
ticles collide, their tangential velocities arc unchanged, 
while the relative normal velocity is decreased by a fac- 
tor of — r, where the negative sign merely indicates that 



they move apart after a collision. Using momentum con- 
servation one can write the final velocities (indicated by 
primes) in terms of the initial velocities 
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In the above equation, the subscript n denotes the ve- 
locity component along the line connecting the centers 
of particles 1 and 2. The energy lost in each collision is 
therefore 
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In this study, we focus on the quasi-elastic limit, i.e. 
1 — r <C 1. Physically, this limit is relevant to hard 
particles such as glass or steel beads. 

Using the above expression for the energy dissipated 
in each collision, one can estimate the sink term, /, the 
mean energy lost per unit area per unit time. Consider 
a particle moving with speed v. During each collision, 
it loses, on average, energy proportional to (1 — r 2 )v 2 . 
In unit time, it collides roughly v/l times. In unit area, 
there are p particles, and consequently 

Ict(l-r 2 )v 2 ^pct(l-r 2 )pT 3 / 2 /l. 



B. The Coefficient of Thermal Diffusivity 

Suppose there is a temperature gradient in the re- 
direction. As a result, there will be an energy flux along 
this axis. To calculate this flux, let us consider the num- 
ber of particles crossing a line perpendicular to this direc- 
tion in a time interval At. We define "crossing the line" 
as having any part of the particle over the line during 
this time interval. To cross the line within At, a parti- 
cle with speed v must have its rightmost point vAt or 
closer to the line. Thus, only grains in an area (d + vAt) 
to the left of the line can pass the line from the left in 
At. In fact, only one half of these particles are mov- 
ing to the right, so the number of particles crossing the 
line per unit cross-sectional length in At is ^p(d + vAt). 
In a steady state, the number of crossing events from 
each side must balance. Any energy flux is due to the 
fact that particles coming from the right are at a dif- 
ferent temperature (Tr) than those from the left (Tj,) 
and thus, qAt = \p(d + vAt)(T L - T R ). Consider two 
grains on opposite sides of the line. Their centers are 
approximately a distance (I + d) apart, so the tempera- 
ture difference is roughly (I + d)dT/dx. The coefficient 
of thermal conductivity is therefore k = —q/(dT/dx) = 
jp(l+d){d+vAt) I 'At. A natural choice for At is the typi- 
cal collision time I /v. This choice is small enough to avoid 
multicounting and is sufficiently large to ensure that heat 
transfer does occur. Our heuristic picture therefore esti- 
mates the thermal diffusivity by 
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This is a rough approximation; the actual prefactors de- 
pend on the velocity distribution of the grains. Therefore, 
we generally assume 



pjal + dfVT 
I 



The energy balance equation (g) takes the form 



d 
dx 



(al + d) 2 ^dT 
I ' dx 



(1 



(4) 



where 7 is the ratio of prefactors in the expressions for 
/ and k. In the following subsections, we discuss how 
to obtain the dimensionless coefficients a and 7 self- 
consistently. To proceed, it is necessary to relate p and 
T through the equation of state. Additionally, the mean 
free path I must be expressed in terms of p. 



Indeed, in the limits of high and low density this expres- 
sion reduces to (|J) and (0), respectively. It is useful to 
compare this expression with the van der Waal's equa- 
tion of state which takes into account long range attrac- 
tion and hard core repulsion [^j. For an inelastic gas 
of hard spheres, there are no long range forces, and the 
van der Waal's pressure for a two-dimensional gas reads 
Pvdw = PPcT/(p c — 2p) for p <C p c - In the low density 
limit the pressure given by Eq. (||) agrees with the van 
der Waal's expression to second order in pj p c . Further- 
more, one can also compare (||) to Tonks' series expan- 
sion |j| for the pressure of a two-dimensional gas of hard 
spheres, P Tonks = pT[l + 1.814(p/p c ) + 2.573(p/p c ) 2 ]/[l- 
1.307(p/p c ) 3 + 0. 307 (p/pc) 4 }, valid for all densities. Over 
the entire density range, the two expressions differ by 
less than 1.3%. In contrast, we found that the van der 
Waal's formula is inadequate for describing the high den- 
sity limit. Hence, we use the interpolation formula (||) for 
the pressure. 



C. The Equation of State 

For the system to be in a steady state, the pressure 
must be constant throughout. The equation of state re- 
lates the pressure and the temperature to the density. 
For example, in the low density limit, the ideal gas law 
holds, PV = Nk B T. Using k B = 1 and p = N/V, one 
has 



P = pT. 



(5) 



On the other hand, in the high density limit, the mean 
free path, which is simply the interparticle spacing, is 
much less than the particle diameter (Z <C d). Denoting 
the close-packing density by p c , one finds 



P_ 
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In this limit, the center of a grain is confined to an area 
of the order of I 2 , so the entropy per particle, S, equals 
ln(/ 2 ) plus a function of temperature. From Eq. (|6|), S 
depends on density only through the term 21n(p c — p). 
Using Maxwell's relation, 
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we obtain the pressure in the limit p — > p c , 

2p 2 T 



P 



Pc~ P 



(J) 



We therefore propose the following interpolation formula 
for the pressure 



P = pT- 



(8) 



D. The Mean Free Path 

The mean free path can be expressed in terms of the 
density and the diameter. In the low density limit one 
has 



I = 



1 



V8pd' 

while in the high density limit, Eq 

, Pc - P 



gives 



2 Pc 



d. 



Again, we use these high and low density limits to inter- 
polate a general expression for the mean free path. Using 
the 2D close-packing value p c = 2/^/~3d 2 , we find 



I 



1 Pc - P 

VSpd Pc-ap' 



(9) 



where a = 1 



E. The density equation 

Eqs. (^) and (^) express the temperature and the mean 
free path in terms of the density. Substituting these 
expressions into Eq. (|4|) yields a second order differen- 
tial equation for p. Using for convenience the variable 
z = p c j p, we have 



dx 



(z 2 + 22I-1) (az(z - 1) + 



f (z- a 



dz 



(z-a)(z-iy/ 2 z 3 / 2 (z + l) 5 / 2 
32 1 - r 2 z-a ~ 



dx 



3d 2 7 + l) 3 / 2 V Z 
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(10) 
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This equation is complemented by the boundary condi- 
tions 



dz 
dx 



and 



x—L 



W / — dx = N, 



(11) 



where W is the width of the system and N is the to- 
tal number of particles. The latter condition merely re- 
flects conservation of particles. The former condition is 
a consequence of the fact that the temperature gradient 
vanishes at the elastic wall. 

In principle, 7 and a are two dimensionless factors 
which can be calculated exactly from the velocity distri- 
bution. However, this distribution is poorly understood. 
Nevertheless, it is still possible to estimate these pref ac- 
tors by comparing the theoretical predictions of Eq. ( |ll| ) 
with numerical simulations in the limit of high and low 
densities. 

In both of these extreme cases the governing equa- 



tion (10) can be solved analytically. In the high den sity 
limit, (p c — p) — > 0, or equivalently, z — 1 <C 1, Eq. ( |1C| ) 
reduces to -£?Vz 



- 1 « 1, Eq. 
It is convenient to 



0, or equivalently, z 

= ^F^ / ^ T ■ 
write this equation in terms of the temperature, T(x) 
From (H), T{x) cx z — 1 when z-l« 1, and the temper- 
ature obeys 



dx 2 



with £ = d\J~ij2(\ — r 2 ). Solving this equation subject 
to the boundary condition dT/dx — at x — L gives the 
temperature profile 



T(x) = T(L) cosh 2 



L — x 



Far from the elastic wall, L — x £, the temperature de- 
cays exponentially in agreement with Haff's calculation 



T(x) = T{0) e - x/i . 



(12) 



wall is implemented in such a way that for r = 1, parti- 
cles have a Boltzmann velocity distribution with an av- 
erage temperature equal to one. Specifically, any par- 
ticle that collides with the wall at x — is ejected 
with a positive v x drawn from the probability distri- 
bution P{v x ) = v x cxp(— v 2 /2) and a v y drawn from 
P(v y ) — exp(— v 2 /2)/\/2tt ||. However, the behavior 
of particles in the bulk of the system is independent of 
the details of the boundary condition (see Fig. || below). 
This is not surprising, but rather a necessity for a ther- 
modynamic theory to be valid. 

Numerical simulations confirm the exponential decay 
in the quasi-elastic limit. Furthermore, the decay length 
£ can be measured for various degrees of inelasticity r (see 
Fig. ||). We verified that indeed £ oc (1 — r 2 ) -1 / 2 as sug- 

2.26 



gested by Eq. (|12|). The value of the prefactor 7 
can be found from the simulations as well. 




FIG. 2. The behavior of the decay length £ in units of par- 
ticle diameters as a function of 1 — r 2 in the high density limit. 
The solid line is 1.06(1 — r 2 )" 1 ^ 2 , and the diamonds are values 
obtained from simulations with an average normalized density 
(p/ Pc) of approximately 0.8. 



We now turn to the low density limit, 
and Eq. @ reduces to £jz 3 / 2 = 
result, 



Here, z 



» 1, 
As a 



Eq. (g) shows that in the high density limit, I cx z—1 cx T. 
Both the temperature and the mean free path decay ex- 
ponentially with the distance from the heated wall. The 
decay length £ is much larger than the mean free path. 
For the continuum description to be valid, £ must also 
be much larger than the diameter of a particle, d, i.e., 



(13) 



Since the prefactor 7 is of order unity, we learn that the 
continuum theory is valid only in the quasi-elastic limit, 
vT — r <C 1. (Note that the restriction \J\ — r <C 1 is 
stricter than 1 — r <C 1.) 

To test the theoretical predictions, we performed an 
event-driven simulation of the system |1Q|. The heated 



dz 
dx 



z - z L 



(14) 



where 77 = 0^/37/64(1 — r 2 ), and we have used the 



boundary condition of Eq. (11) and the notation Zl 
z(L). Simulations show that the density near the heated 
wall is significantly smaller than the density at the elastic 
wall. Thus, for x <C L and z > zj,, we find the following 
approximate solution 



Z{x) W Zq 



V 



The inverse of the density, and hence the temperature, 
decays linearly with x in the low density regime near 
the heated wall. That the low density decay length 
77 < — > (1 — r 2 ) -1 ' 2 is similar to the high density decay 
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length £ ~ (1 — r 2 ) -1 / 2 , reflects the fact that the under- 
lying differential equation is second order. 

In the nearly elastic case (1 — r <C 1), zl is close to z, 
and Eq. (|lj) can be integrated exactly 

- = - \\fz\fz - z L + z L ln(</z + \/z- z£)] z } x) . 

Comparing this prediction with low density numerical 
simulations allows us to determine r\ at various r values 
(see Fig. ^ and verify that indeed 77 oc (I — r 2 ) -1 / 2 . The 
constant of proportionality and the previously calculated 
7 yields a = 1.67. 




FIG. 3. The behavior of the decay length 77 in units of par- 
ticle diameters as a function of 1 — r 2 in the low density limit. 
The solid line is 0.542(1 — r 2 ) _1//2 , and the diamonds are val- 
ues obtained from simulations with an average normalized 
density (p/pc) of approximately 0.1. 



In the low density limit, the equation of state P = pT, 
implies that T(x) — —z(x). For the state to be locally 
very close to equilibrium, the mean free path must be 
much less than the length scale over which the tempera- 
ture is changing, i.e., I <C r\p c jP. This condition is 



V / l-^ 2 «a^, (15) 

when the temperature is of order unity. Since a and 7 
are also of order unity, the hydrodynamic description is 
again valid only in the quasi-elastic regime (VI — r <C 1). 

Our treatment so far has concentrated on cither the 
high or the low density limit, where analytical expressions 
were possible. For systems that include both high and 
low density regions, Eq. (jic]) can be solved numerically 
using the boundary and normalization conditions (|ll| ) 
and the previously calculated values for 7 and a. By ex- 
amining a number of simulations that include a range of 
densities we determined that the optimal values for the 
prefactors are 7 = 2.26 and a — 1.15. While this value of 
7 is consistent with the value obtained in the high den- 
sity calculation, this a value is slightly lower than our 
prediction. A typical system with r = 0.99 is shown in 
Fig. ||, and it is seen that the predictions of hydrodynamic 
theory match the numerical data over wide density va ri- 
ations. For this simulation, the ratio \/2(l — r 2 )/~/a 2 is 



approximately 0.1, which is not of order unity, so the 
condition of Eq. (jig) is satisfied. 




50 100 150 200 250 



position (x/d) 

FIG. 4. Comparison of a simulation (plus signs) of 1000 
particles at r = 0.99 with the numerical solution (solid lines) 
to the full differential equation ( |lo| ) with 7 = 2.26 and 
a — 1.15. The temperature is not one at the heated wall 
due to the effect discussed in Section III. The horizontal axis 
is the distance from the heated wall in units of the particle 
diameter (d). 



III. NONEQUILIBRIUM STEADY STATE 



The calculation in the previous section assumes that 
the steady state is very close to equilibrium, and that 
pressure and temperature can be used to describe the 
system. This requires that condition ( |i~3| ) is satisfied in 
the high density limit, or condition ([l5]) is satisfied in the 
low density limit. 

An important question is: when does the behavior 
of the system changes qualitatively? One such transi- 
tion occurs when r becomes low enough for the system 
to undergo inelastic collapse. Here strong correlations 
and large density variations develop, and applying hy- 
drodynamics becomes impossible. Even when r is slightly 
higher than the critical value for inelastic collapse, the at- 
tractors mentioned in the introduction may still be strong 
enough to build correlations. Where this breakdown oc- 
curs is determined by the degree of inelasticity, the den- 
sity and the total number of particles Jll| . For each pair 
of values of density and total number of particles, there is 
a value of r which divides two different kinds of behavior: 
loose sand and coherent sand. 

However, the theory developed in the previous section 
is strictly for the quasi-elastic limit, 1 — r <C 1, so there 
may exist systems that, although elastic enough to avoid 
inelastic collapse, still have r far enough from one that 
the hydrodynamics do not apply. Since correlations be- 
tween particles are built up through inelastic collisions, 
high density regions are more liable to inelastic collapse 
prf, while for low density regions, correlations are harder 
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to establish. Therefore, we will investigate the low den- 
sity limit in order to observe the breakdown of the hy- 
drodynamic description as the degree of inelasticity in- 
creases. 

When the temperature variation within a mean free 
path is significant, the system is unable to reach local 
equilibrium. Therefore, particles carrying energy away 
from the heated wall cannot share this energy with the 
slow particles returning from the higher density region 
near the elastic wall. This inefficient mixing leads to a 
temperature gap - the average energy of particles with 
v x > is greater than that of the particles with v x < 
(see Fig. ||). Furthermore, near the heated wall, the tem- 
perature drops by approximately 15% over a mean free 
path, which suggests that such a system will be unable 
to reach local thermal equilibrium (see Fig. ||). 




FIG. 5. The difference in temperature for particles with 
positive velocities (solid line) and with negative velocities 
(dashed line). The data is from a simulation of 1500 particles 
with r = 0.95 and total area fraction 0.01. 




FIG. 6. The mean free path, i.e., the average distance until 
the next collision for particles at a given position. The dip 
near the heated wall is due to collisions with the wall. The 
data is from the same simulation used in Fig. ^. Both the 
mean free path (I) and the position (x) are in units of particle 
diameters (d). 

The probability distribution functions for v x (the ve- 



locity perpendicular to the heated wall) provide an illu- 
minating way to measure this deviation from equilibrium 
(see Fig. ^). Note the asymmetry: the v x > tail is 
longer than the v x < tail. This is consistent with our 
understanding that the v x > particles have more en- 
ergy than the v x < ones. A rough quantification of the 
deviation from equilibrium is provided by condition (flF 



when \/2(l — r 2 )/ja 2 is of order unity, the theory breaks 
down. For r = 0.95 (the value used in the Figures), this 
quantity is approximately 0.3. Note that this condition 
involves r only and is not sensitive to the local density 
p. This suggests that the behavior of regions of the sys- 
tem with different densities should be similar. In fact, 
we observed that the velocity distribution function obeys 
scaling (see Fig. H), i.e., 



P(v x ,x) 



1 



g(x) \g(x) 



(16) 



The function <fr is independent of boundary conditions 
(see Fig. ||). Thus (v x (x)} cx g n (x), where the constant 
of proportionality depends only on the shape of the func- 
tion 4>. Note that this shape will depend on r, the degree 
of inelasticity; the more inelastic the particles are, the 
more skewed the velocity distribution is. The probability 
distribution for the y-componcnt of the velocity behaves 
similarly: 



P(Vy,X) 



1 



g{x) \g{x) 



(17) 



but here ip is a symmetric, nearly Gaussian, function. 
Additionally, the same velocity scale g(x) characterizes 
the transverse and the longitudinal velocity distributions. 
Therefore, while at each position x there is no longer a 
single hydrodynamic temperature, there is a well defined 
characteristic velocity scale, g(x), so that the granular 
temperature, \(v 2 x + Vy), is proportional to g 2 (x). This 
scaling suggests that, although the system may deviate 
significantly from equilibrium, it can still be treated using 
some of the tools of the previous sections. 

Specifically, return to the equation for energy bal- 
ance (0): dq/dx = —I, where / is the energy loss due 
to collisions per unit time per unit area, and q is the en- 
ergy flux, the heat transfer per unit cross-sectional length 
per unit time. As discussed previously, 



/cx (1 



r 2 )dp 2 v 3 



oc p 2 g 3 , 



while the heat flux is approximately 

1 = {{\pv 2 )v x ) oc pg 3 . 

Conservation of momentum flux suggests that pg 2 is con- 
stant, and the energy balance equation gives 



dg 
ax 



(18) 



G 



which indicates that g 2 , and hence the temperature, de- 
pends linearly on x. This prediction is consistent with 
our simulational results (Fig. W. 
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FIG. 7. The probability distributions for v x at various posi- 
tions x = L/5, L/2, 4L/5, where L is the length of the system 
(1330 particle diameters in this case). The simulations are for 
a system of 1500 particles with area fraction 0.01 at r = 0.95. 
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FIG. 8. The scaling distribution 4>(z) — P(v x ,x)g(x). 
The data of Fig. W have been rescaled according to 
Eq. (^), so three sets of data are from various positions, 
x — L/5(+), L/2(*), 4L/5(o), in a system with the usual en- 
ergy enput at x — 0. The plot also includes probability dis- 
tributions from a simulation of a system that is identical to 
that used in Fig. ^ but with a different type of forcing at 
the heated wall. Data from this simulation is displayed for 
x = 7L/10(A), 4L/5(n). 



IV. CONCLUSION 



no longer close to equilibrium in the low density limit. 
However, the scaling behavior of the velocity distribu- 
tions suggests that a hydrodynamic treatment can still 
be useful in describing the system. 

In this non-equilibrium regime, we found that the par- 
ticle velocity distributions were non-Gaussian. Indeed, 
deviations from normal distributions have been observed 
in theoretical Jl^,|l2|,^3) and experimental studies [fl4|| . 
In addition to this variation in the velocity distributions, 
systems of this sort can produce highly inhomogeneous 
spatial distributions, as has been noted elsewhere in one 
and two j?],|l5| dimensions. A recent experimental 
study examined the spatial distribution of hard particles 
in two dimensions in the presence of an energy input |l6| . 
Their data is in qualitative agreement with our theoreti- 
cal predictions, and inhomogeneous spatial distributions 
reminiscent of Fig. (1) are observed. 
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FIG. 9. Characteristic velocity scale as a function of posi- 
tion. Here we plot {(v x + Vy) 1 ^ 2 ) oc g(x) (the dashed line). 
This data was obtained from the simulation used to make 
Fig. [7 and the values for g(x) were used to do the rescaling in 



Fig. 



We also plot («| +Vy) oc g 2 (x) (the solid line) to show 



that it is a roughly linear function of position, as predicted 
by Eq. (|lS 



In this work, we examined the steady-state behavior of 
a weakly inelastic two-dimensional driven granular sys- 
tem. We found that a hydrodynamic formulation pro- 
vides a satisfactory description of the near-equilibrium 
behavior of the system in the quasi-elastic limit. We ex- 
tended Haff 's theory to the low density limit and found 
that the corresponding temperature profile varies linearly 
in space. For slightly higher inelasticities, the system is 



[1] H. M. Jaeger, S. R. Nagel, and R. P. Behringer, Phys. 

Today 49, 32 (1996). 
[2] P. K. Haff, J. Fluid Mech. 134, 401 (1983). 



7 



[3] J. T. Jenkins and M. W. Richman, J. Fluid Mech. 192, 
313 (1988). 

[4] Y. Du, H. Li, and L. P. Kadanoff, Phys. Rev. Lett. 74, 
1268 (1995). 

[5] E. L. Grossman and B. Roman, Phys. Fluids 8, 3218 
(1996). 

[6] J. Maddox, Nature 374, 11 (1995). 

[7] S. E. Esipov and T. Poschel, "Boltzmann equation and 

granular hydrodynamics," preprint. 
[8] H. Hayakawa, S. Yue, and D. C. Hong, Phys. Rev. Lett. 

75, 2328 (1995). 
[9] L. Tonks, Phys. Rev. 50, 955 (1936). 
[10] I. Goldhirsch and G. Zanetti, Phys. Rev. Lett. 70, 1619 



(1993). 

[11] S. McNamara and W. R. Young, Phys. Rev. E 53, 5089 
(1996). 

[12] Y. D. Lan and A. D. Rosato, Phys. Fluids 7, 1818 (1995). 
[13] Y. H. Taguchi and H. Takayasu, Europhys. Lett. 30, 499 
(1995). 

[14] S. Warr, G. T. H. Jaques, and J. M. Huntley, Powder 
Tech. 81, 41 (1994). 

[15] S. McNamara and J.-L. Barrat, "The energy flux into a 
fluidized granular medium at a vibrating wall," preprint. 

[16] A. Kudrolli, M. Wolpert, and J. P. Gollub, "Cluster for- 
mation due to collisions in granular material," preprint. 



8 



